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Abstract 

In 1930 Sergey L. Sobolev 0S has proposed a construction of the solution of 
the Cauchy problem for the hyperbolic equation of the second order with variable 
coefficients in 3-d. Although Sobolev did not construct the fundamental solution, 
his construction was modified later by Romanov 0! to obtain the fundamental 
solution. However, these works impose a restrictive assumption of the regularity 
of geodesic lines in a large domain. In addition, it is unclear how to realize those 
methods numerically. In this paper a simple construction of a function, which is 
associated in a clear way with the fundamental solution of the acoustic equation 
with the variable speed in 3-d, is proposed. Conditions on geodesic lines are not 
imposed. An important feature of this construction is that it lends itself to effective 
computations. 
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1 Introduction 

In 1930 Sergey L. Sobolev, one of the most distinguished mathematicians of the 20 st 
century, has published two papers ms where he has constructed the solution of the 
hyperbolic equation of the second order in the 3-d case with variable coefficients in the 
principle part of the operator. It was assumed that these coefficient depend on spatial 
variables. This result is readily available in the textbook of Smirnov [6]. Sobolev did 
not find the fundamental solution. The main reason of this was that the notion of the 
fundamental solution was unknown in 1930. Still, in his books about inverse problems 
00 Romanov has modified the method of Sobolev to construct the fundamental solution 
of that equation. However, the constructions of both Sobolev and Romanov impose a quite 
restrictive assumption on the variable coefficients in the principal part of the hyperbolic 
operator. Let the time variable t e (0, T). It is assumed in above cited publications that 
geodesic lines generated by these coefficient are regular in a domain Q (T) C R 3 . The 
larger T is, the larger Q (T) is. So, Q (oo) = R 3 . The regularity of geodesic lines in the 
domain Q (T) means that for any two points x, y G Q (T) there exists a single geodesic 
line connecting them. 
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Another construction of the fundamental solution of that equation can be found in 
the book of Vainberg pi. The construction of pi imposes the non-trapping condition 
on variable coefficients in the principal part of the hyperbolic operator. In addition, the 
technique of [10] relies on the canonical Maslov operator, which is not easy to obtain 
explicitly. Thus, even though the structure of the fundamental solution in any of above 
constructions can be seen, formally at least, still many elements of this structure cannot be 
expressed via explicit formulas. For example, neither the solution of the eikonal equation 
of the method of Sobolev, nor the Maslov construction cannot be expressed via explicit 
formulas. 

The numerical factor is important nowadays. However, it is not immediately clear how 
to compute numerically fundamental solutions obtained in above references. Therefore, it 
is also unclear how to compute solutions of Cauchy problems for heterogeneous hyperbolic 
equations if using those fundamental solutions. 

In this paper we propose a simple method of the construction of a function, which is 
associated in a clear way with the fundamental solution for the acoustic equation in the 
3-d case with the variable coefficient. We impose almost minimal assumptions on this 
coefficient. Geodesic lines are not used. It is important that our function can be both 
accurately and effectively approximated numerically via the Galerkin method as well as via 
the Finite Difference Method. Thus, we show that this function can be straightforwardly 
used for computations of the Cauchy problem for the heterogeneous acoustic equation. 

The author was prompted to work on this paper while he was working on publications 
BE] about reconstruction procedures for phaseless inverse scattering problems. Indeed, 
in BE] the structure of the fundamental solution of the acoustic equation in time domain 
with the variable coefficient in its principal part is substantially used. 

2 Construction 

Below x G R 3 and ff C R 3 is a bounded domain. Let cq and c\ be two constants such 
that 0 < c 0 < Ci. We assume that the function c(x) satisfies the following conditions 


c E C 1 (R 3 ) , c G [co, Ci] , (2-1) 

c (x) = 1 for x G R 3 \f2. (2.2) 

Consider the following Cauchy problem 

c (x) ytt - A x y = f (x, t ), (x, t)Gl 3 x [0, oo), (2.3) 

y |t=o= Vt |t=0= o. (2.4) 


It is assumed here that / is an appropriate function such that there is a guarantee of 
the existence of the unique solution y G H 2 (R 3 x (0,T)) ,VT > 0 of this problem, see, 
e.g. Theorem 4.1 in §4 of Chapter 4 of the book of Ladyzhenskaya [3] for sufficient 
conditions for the latter. We will specify / later. The fundamental solution for the 
operator c (x) d f — A is such a function P (x, £, t, r) that the solution of the problem (12.3)) . 
(12.4)) can be represented in the form 

OO 

y(x,t) = J j P(x,£,t,r) f (£,T)d£dT. (2.5) 

0 R3 
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Below we modify formula (12.51) . 

Let £ € K 3 and r > 0 be parameters. Consider the following Cauchy problem 

c(x) uu = A x u + S(x - f) 5 (t - t) , (2.6) 

u |t=o= Ut |t=o— 0. (2-7) 

2.1 Heuristic part of the construction 

It is convenient for us to work in this subsection with a purely heuristic derivation. Rep¬ 
resent the solution of the problem (12.611 . (12.71) as u — u 0 + v, where u 0 is the fundamental 
solution of the wave equation, 


u 0 = 


5(t - t - \x - g|) 
4-7T \x — £| 


Thus, 


d^u 0 = A x u 0 + 5(x-£)5(t-T), 


Uq |t=o— Uot |t=o— 0. 

Hence, the function v satisfies the following conditions 

|x-£|) 


c (x) v tt = A x u - (c (x) - 1) 


4-7T \x — £| 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


( 2 . 11 ) 


v |t=o— Vt |t=o— 0. (2.12) 

Consider the operator A, 

t 

Mf) = Jf(y)dv 

0 

for appropriate functions /. Purely heuristically again apply the operator A 4 to both 
sides of equation (12.lip . Denote w (x, £, t, r) = A 4 ( v ). Then (12.111) and (12.12(1 imply that 


c (x) w u - A x w = - (c (x) - 1) ^ r . ^ ^ H (t~- t - \x - £|), 

4-7T \x — £| 

U! |t=o= u>t U=0= 0, 

where H ( 2 ) is the Heavyside function, 


H(z) 


1, z > 0, 

0, z < 0. 


(2.13) 

(2.14) 
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2.2 Rigorous part of the construction 

Starting from this point, we are not acting henristically anymore. To the contrary, we 
work rigorously everywhere below. 

Consider the Cauchy problem (12 .1 3f) . (12. 14)1 . Denote 

9 {x, f, t,r) — — (c (x) - 1) ^ T . ^ ^ H (t - t - \x - f |). 

47T \x — £| 

It follows from (12.2)1 that for any hxed pair (£, r) G R 3 x [0, oo) and for any T > 0 
functions g,gt G L 2 (R 3 x (0,T)). Hence, Theorem 4.1 of §4 of Chapter 4 of the book 
of Ladyzhenskaya m, m, o as well as other results of that chapter imply that 
for any hxed pair (£, r) G R 3 x [0,oo) and for any T > 0 there exists unique solution 
ui G H 2 (R 3 x (0, T)) of the Cauchy problem (j2. 1311 . (12.141) . Furthermore, it was shown in 
the proof of that theorem of [3] that this solution can be effectively constructed numerically 
via the Galerkin method. It is also well known that it can be numerically constructed via 
the Finite Difference Method. In addition, the energy estimate implies that 

w (, x , £, t, r) — 0 for t < r. (2.15) 

Also, it follows from (12.13)1 that w — w (x,£,t — r) 

Consider now the function w defined as 

w (x,£,t — t) = w (x, £,t-T) + A 4 (u 0 ) . 


Hence, 


w (x, f, t - r) = w (x, £,t-r) + — T j X ^ H (t - r - \x - ^|). (2.16) 

6 • 47T \x — £| 

Using (l2.8D - fl2.10ll . (12.13)1 . (12.14)1 and (j2. 16)1 and applying direct calculations, we obtain 
that the function w satisfies the following conditions 

c (x) w tt - A x w = 6(x-£) — ^ H (t - r) , (2.17) 
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w |t=o= w t |i=o= 0. (2.18) 

Consider now an arbitrary function / (x, t) such that 


/ (x, t) G C 4 (R 3 x [0, 00 )) , 


df f (x, 0) = 0, k = 0, 1, 2, 3,4, 

/ (x, t ) = 0, Vx G R 3 \G/, 

where Gf is a bounded domain depending on the function /. Consider the function 
Pf (x, t ) defined as 


Pf ( x, t ) 



0 G f 


T) f (€, T ) d i dT - 


(2.19) 
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Using the integration by parts, (I2.15p and (I2.16p . we obtain that there exist four derivatives 
of the function pj ( x,t ) with respect to t. In particular, 

t 

dtPf OM) = J J w (u €,t-r) d 4 f (f, r) d£dr. (2.20) 

0 G f 

Next, applying the operator c (x) d 2 — A to both sides of (I2.20p and using (I2.15P and 
(12. 1711 . we obtain 


(c (x) d 2 - A) dfp f (x, t) = j &—SLg}f ( Xj r ) dr. (2.21) 

o 

Using integration by parts in (I2.2ip as well as (??), we obtain 

(c (x) d 2 t - A) d?p f (x, t) = f (x, t ). (2.22) 

In addition, it follows from (12.151) . (I2.16K and (12.20p that 


d?Pf U=o= dtPf U=o= 0. 

Since the function d 4 pf G H 2 (R 3 x (0, T)) , VT > 0, then the uniqueness theorem for the 
problem (12.31) . (12.41) and (12.201) imply that 


y(x,t) 


J J w(x,€,t-T)d*f (£,t) d^dr. 

0 G f 


(2.23) 


Thus, the solution of the problem (12.31) . (12.4p is constructed in (12.231) using the function 
w. The formula (12.231) is our first analog of the formula (12.51) . The integration by parts 
in (I2.23P and the use of (I2.15P and (I2.16P leads to the second analog of the formula (12.51) . 


y(x,t ) 


J J d 2 w(x,C,t 

0 G f 


r)d 2 f (£ 7 r)d£dT. 


(2.24) 


2.3 Numerical comments 

Recall that the function w can be accurately numerically approximated via either the 
Galerkin method or the Finite Difference Method. The derivatives of the function / 
can be found analytically, if / is given by an explicit formula. However, if / is given 
with a noise, then a regularization method should be applied, see, e.g. the book of 
Tikhonov and Arsenin [Qj. In particular, it is explained in this book how to stably 
differentiate noisy functions using the regularization. Furthermore, numerical examples of 
stable computations of first and second derivatives are presented in [9]. The form (I2.24p 
might be sometimes more convenient than the form (12.231) since second derivatives of 
noisy functions are obviously more stable to calculate than fourth derivatives. 

In summary, formulas (12.231) and (I2.24P imply that the solution y (x, t) of the problem 
(12.31) . (12.4p can be effectively numerically calculated via three steps: 
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1. Step 1. Compute the function w either via the Galerkin method or via the Finite 
Difference method. 

2. Step 2. Compute corresponding t— derivatives of the function / should be computed 
either analytically, if / is given by an explicit formula, or numerically using the 
regularization, if / is given with a noise. 

3. Step 3. Apply one of formulas (I2.23j) or 02.241) . 
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